Implementation of quantum operations on single photon qudits 
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We show that a general linear transformation from one single photon qudit to another, the di- 
mension of which can be either equal or unequal to that of the first one, can be implemented by 
linear optics. As an application of the scheme we elaborate a method to deterministically realize 
any finite-element Positive Operator Value Measure (POVM) on single photon signals, which is also 
generalizable to any quantum system in principle. 
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Linear optics is considered as one of the promising can- 
didates for quantum computing (for a recent overview 
see, e.g., [l|) and can be also applied to many other ar- 
eas such as quantum cryptography (for a review see. e.g., 
In these applications an essential technique is the 
implementation of all possible operations, including gen- 
eralized quantum measurements in the form of Positive 
Operator Value Measures (POVMs), on the signals en- 
coded as photon states by practical linear optics circuits. 

A typical and important case of the signal states is 
single photon qudits, i. e., the linear combinations of the 
modes a|.|0), k = 1,---,N (multiple-rail encoding). It 
was proposed by Reck et al. 3] that any unitary operator 
U G U{N) on the A^-dimensional qudits, Cia||0), 
can be realized by an N-port interferometer, which is 
an array of beam splitters and phase shifters performing 
SU{2) elements, because this unitary operator can be 
decomposed into the product of these SU{2) elements 
(see Fig. 1). This scheme was further studied in @| 
and has been applied to a variety of research fields in 
quantum information theory and experiment. 

The generalization of the scheme is the implementa- 
tion of all possible linear maps on single photon qudits, 
which is a fundamental task in processing quantum in- 
formation. It is intimately related to the realization of 
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FIG. 1: The unitary transformation module constructed 
with beam splitters (dark square) and phase shifters (white 
square). Any unitary operator U can be decomposed into the 
product of SU (2) elements implemented by the beam splitters 
and the phase shifters, and the maximum number of beam 
splitters needed is N{N — l)/2. The input ports are with 
unprimed numbers and output ports with primed numbers. 



POVMs that are at the heart of many quantum infor- 
mation processing protocols. A finite-element POVM is 
a set of non- negative operators {11^}, where 11^ are its 
elements, satisfing 
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with / being the identity operator. It has been proved 
that any rank-one POVM in the form of 11,; = kf\(l)i){<j)i\, 



where 



i) ^ij ^'^d \ki\ < 1, can be realized by the 



Neumark extension 7], which extends the POVM ele- 
ments to the orthogonal projectors in a larger space. For 
the input signals prepared with single photons, such a 
POVM can be implemented with linear optics circuits, 
performing unitary transformations, and photon detec- 
tors only [81]. The realization of POVMs with arbitrary 
rank is, however, much more difficult. Since Hi > 0, it 
can be decomposed into Hi = AjAi with Ai being upper 
triangle matrix The general POVM will be imple- 
mented if we simultaneously realize the maps. 
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and successfully detect these outputs. The detection op- 
erators Ai (and their transforms by an arbitrary unitary 
operator A[ = UiAi) can be any allowed linear map in 
quantum mechanics with equal dimensional input and 
output signal space, i.e., an N x N square matrix. 

The detection operators of a POVM are, however, only 
a special case of more general maps called quantum op- 
erations (QOs). A QO connects pair of input and output 
states via the map. 
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£ is a linear, trace-decreasing map that preserves the 
complete positivity (CP), and generally occurs with non- 
unit probability Tr(£(pi„)) < 1. The general form of £ 
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IS given as 
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with the Kraus operators Ki satisfying the bound 
'Et KjKi < I. If a QO transforms pure states to pure 
it is called a pure map. In this case there is only one 
term K pinK^ in the above equation with K being a con- 
traction, i.e., \\K\\ < 1. For a single pure state input 



Pin = |i/'m)(V'm| [III j the QOs in Eq. (3) can be there- 
fore written in the form 11211. 
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The output signal of such a linear map, with K as con- 
traction, is detected with a probability {ipin\K^ K\ipin) < 
1. These contractions can be more general than the de- 
tection operators of a POVM because the input space 
dimension A^i and the output space dimension N2 of K 
can be different, so K is an A''2 x Ni matrix whose entries 
are complex numbers. The detection operators Ai of a 
POVM correspond to a special case of contractions when 
Ni^ N2= N. 

We now address the problem of how to realize any 
possible linear map K on single photon qudits with only 
three unitary operator modules of the kind shown in Fig. 
1. We realize K by its unitary dilation, U, the unitary 
operator constructed from K in a larger space, which we 
obtain by using the direct sum extension of the system 
with an ancilla, Tis © Ha- In terms of Hilbert space 
dimensionality, this scheme minimizes the physical re- 
sources needed to realize a QO [l3| . We embed the state 
vector (ci, C2, • • • , catJ-^ (T stands for transpose) of the 
input signal, Itpin) = '^j'^i ^^^^ ^ larger space 

and map it by W to a vector containing the state vector 
of the output, iV'out) = '^f^i CiO.^^ \0) (unnormalized) , of 
K: 
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It should be noted that we realize this unitary dilation 
always with a vacum state ancilla (0, 0, • • • , 0)-^ added to 
the input state vector as a direct sum. Next, we will 
derive the algorithm to generate the unitary dilation U 
of linear map K. 

K^K and KK'' are positive matrices of A^i x A^i and 
N2 x N2, respectively. Suppose A^i > N2, we choose to 
diagonalize KK'' by a unitary operator U. We can also 
construct an A'^i x A^i unitary matrix V to obtain the 
Singular Value Decomposition (SVD), K — UI^V'' , with 
the uniquely determined singular values on the diagonal 



of the N2 X Ni matrix S. If Ni < N2, on the other hand, 
we will choose to diagonalize K'^ K to get a similar result. 

Then we extend the rectangular matrix S to the fol- 
lowing max{Ni, N2) X max{Ni, N2) square matrix 
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where the singular values Ui satisfy \ai\ < 1 since is a 
contraction. The extension of E in the case of iVi < N2 
also takes the above form except that IctatsI is replaced 
by \o'Ni\- Using the fact that S' is still a contraction, we 
can obtain a max{2Ni,2N2) x max{2Ni,2N2) unitary 
dilation 
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S' (/ - E'2)l/2 

(/„I]'2)l/2 

of it (see, e.g., Ex. 1.3.6 in [^), which acts on a space 
H(BH with H being max{Ni, N2) dimensional. We also 
extend U and V to max{2Ni,2N2) by max{2Ni,2N2) 
matrices by adding the identity matrix / in the diagonal 
and zero matrices off the diagonal. A general linear map 
K is therefore realized by the following unitary dilation: 



U = UGV''. 



(9) 



In our setup, we perform its equivalence by acting 
on the spatial mode vector {a\ , • • • , ajy^ ) . The cir- 
cuits to implement V and are the corresponding A'^i- 
port and N2-poTt modules. After the input spatial mode 
vector is processed by V, we redirect the output to a 
TOaa;(2Afi, 2iV2)-port module of G with the input ports 
numbered from iVi -f 1 to max{2Ni, 2N2) in Fig. 1 black 
or a vacum state. Here is some detail about the step to 
implement E through its unitary dilation G. Picking out 
the entries containing only one of the singular values at 
from the matrix of G, we form a 2 x 2 submatrix, which 
can be transformed by a rotation 7i j^,„Q2,(jVi,Af2) to a 
diagonal one: 



cos 6i — sin 9i 
sin 9j cos 6i 



(10) 
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With a series of rotations in the form of Tj_^i+max(Ni,N2) ® 
/rest, for i = 1, • • • , min{Ni, N2), G can be realized by 
min{Ni, N2) beam splitters with the reflection coeffi- 
cients R = 1 ~ and phase shifters giving rise to e''^. 
Therefore, the upper bound of the total number of the 
beam splitters required in the scheme is 
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FIG. 2: The circuit to perform the unitary dilation U oi K, 
the general linear transformation on a pure state input pin ~ 
\tpin){tpin\- We wiU obtaiu two outputs, pout of K and an 
extra p^x, from the corresponding terminals. Because the 
linear map K is not generally inversable, we need to add an 
ancilla pex to pout, i.e., p'i„ — pout © Pex, if we are to convert 
it back to the original pi„ by using the same circuit from the 
inverse direction. 



which is determined by the dimensionahties of the in- 
put and output Hilbert spaces. In the whole extended 
space, we will obtain two outputs after the action of 
the three unitary operator modules: one is the ex- 



act output {a-^ 



^N2 



) of the linear map K from the 



output ports of , and the other is an extra output 
iaN2+v---^"-lax(2Nu2N.2)) ^^1° output ports of G 

numbered from {N2 + 1)' to (maa;(27Vi, 27V2))'. Fig. 2 
displays the scheme that realizes the effect of K on the 
input state pm- 

This linear optics scheme can be directly applied to 
where we need non-unitary transformation on photon 
states, e.g., in the production of single photon qudits 



in any form of X^i^iaJlO), where J2i 



< 1 (possi- 



bly unnormalized) , by multiple-rail encoding, and in the 
enhancement of the entanglement of a pair of partially 
entangled photons like Cifl- |0)ia||0)2, with different 
\ci\ unequal, by one party operation. 

Now we look in some detail at the realization of 
POVMs as one important application of our scheme. We 
start with the simplest situation of n = 2, where the two 
POVM elements are always commutative, [111,112] — 0. 
Suppose that the dimension of the signal space is iV, 
and the N x N detection operators Ai of the POVM 
can be factorized by SVD as Ai — ViYiiUi with Ui, Vi 
unitary and diagonal. We first set up a A^-port mod- 
ule for Ui and, after the signal leaving Ui module, we 
process it with a 2A^-port module to implement the uni- 
tary dilation of Si. From its output ports numbered 
from 1' to N' we get an output IV'mid) ~ 5^iC^i|V'm), 
while from the ports numbered from (TV -I- 1)' to 2N' an- 
other output IV'mid) SicC/ilV'm), with T^lfj = / - S^. 
Then we will just redirect them to modules of Vi and V2 
and finally obtain the outputs Ai|?/'i„)/||Ai|^i„)|| and 
^2|V'Mi)/ll^2|V-'iri)|| fi'om the corresponding terminals. 

For a POVM with the number of elements n > 3, the 
situation is much trickier. Instead of 112, what we realize 
from the corresponding output ports of \4'mid) *^P" 
erator /— Hi . By diagonalization all elements of a general 
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FIG. 3: The setup to perform any POVM with three ele- 
ments. The seven unitary operator modules, two of which 
perform the pairs of operators Ei and Eic — {I — Ef)^''^, 
E2 and E2C = ~ ^2^y^^, respectively, are designed with 
the POVM elements. The detectors at the terminals effect a 
dephasing to eliminate the interference between Ai\tpin) [l8l |. 
and capture the outputs pout,i = AiPi„Al/Tr{AipinAl), for 
i = 1,2,3, with the prbabilities pi = Tr{Aipi„Al). The out- 
puts redirected to only one set of detectors form a probability 
distribution as the mixture, J2iPiP°^t-i ^ J2i^'P^"M US' 
there. 



POVM can be factorized into 11, 



- rrtv2 



UjT,fUi, where the 



different Ui do not generally commute, i.e., [Ui,Uj] ^ 
for i ^ j. In the realization of 112, therefore, we need to 
consider two different situations: (1) If ||ni|| < 1 
because / — Hi > 112 when n > 3, we can find a diagonal 
matrix with HSjH < 1 15[ and a unitary operator 
U2L such that 

n2 = UIJ:icU^^^^*2''U2L^1cUi = U^2^lU2. (12) 

These two matrices and U2L are obtained by a stan- 
dard digonalization procedure following the above equa- 
tion. Then, after performing U2L with a TV-port module, 
S2 as a contraction map can be implemented by a 2TV- 
port linear optics module with at most A'^ beam split- 
ters. To realize A2 completely, we add one more module 
of a proper ¥2- (2) If ||ni|| = 1, after the signal goes 
through the part of circuit implementing I — Hi, some of 
the output ports will be black because the cor resp onding 
components have been projected out by Hi [1^. Then 
we will just inverse the remaining (TV — D) x (TV — D) 
non-zero part of Sic matrix, where D is the multiplicity 
of the unit eigenvalue of Hi, in finding U2L and E2 of 
this size in Eq. (12). 

Repeating the above procedure from the output ports 
where the operator 7 — 111—112 is realized, we add all 
the corresponding modules performing UnL, S* , etc., for 
n > 3, to implement the remaining A3, A^, An, 
respectively. The total number of modules of Fig. 1 
needed in our scheme to realize a general POVM with 
n elements is 3n — 2. As an illustration of this general 
scheme. Fig. 3 shows the setup to perform any POVM 
with three elements. 

We have reduced the problem of realizing a POVM 
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to that of finding a sequence of unitary operators with 
the POVM elements and reahzing them with the ancilla 
states of vacum and then detecting the outputs with 
the standard projective measurements on the extended 
space. The algorithm to obtain these unitary operators 
is given, and the implementation of any POVM, with el- 
ements of arbitrary rank, can be therefore realized for 
single photon input signals. Using this method, we will 
directly obtain the output states of a POVM, which can 
be tailored by choosing the appropriate Vi modules, from 
the corresponding terminals where the signal detectors 
are placed. If our signals are just single photon polariza- 
tion qubits, \ipin) = ci\H) + C2\V) (polarization modes 
H and V) , we can use much simpler circuit to implement 
any POVM on them as in [l3|, where a POVM is re- 
alized as the decomposition of an identity operator but 
the necessary algorithm to obtain, e.g., S* , UnL for the 
implementation of all specified 11^ is not given. Given 
the beyond-linear-optics methods to implement unitary 
operations on more complicated quantum systems than 
single photon states, this scheme can be applied to more 
general situations of photonic states as well as the signals 
of any other type of radiation. 

In the whole extended output space Tir, the overall 
output of the unitary map U implemented by the POVM 

circuit is lij'^^t) = {{Al\iJ,r^)V,■■_AAr^\1p,n)V)'^■ Ap- 
pling a proper dephasing map [l8| to this output as 

'^{p'out)i '^ill obtain a direct sum of AiPinAj in Ht, 
which can be transformed by a contraction map £ [Toj 
to W{pin) = (SiLi A^pinAl)/n in one subspace of Ht- 
The output state of the total map W on a pure state in- 
put \ipin), normaHzed as in Eq. (3), is pout = J2t ^iPinA] 
of a QO mapping \Tpin) from a pure to mixed state. A 
general QO can be therefore realized by a corresponding 
combined map W. To a single set of detectors that ef- 
fect a dephasing map by detecting Ailtpm) coming from 
different terminals, what is being measured is effectively 
the output of a general QO (see the caption of Fig. 3). 

In conclusion, we have presented the linear optics 
schemes (including the photon detection) to realize all 
QOs and POVMs on a single photon qudit. The circuits 
to perform all the relevant tasks are only the combina- 
tions of some scalable unitary operator modules which 
have been widely applied in quantum information pro- 
cessing. Given current technologies, our schemes can re- 
alize any permissible map on single photon signals in a 
deterministic way. 
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